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SUNRISE/IMAX OBSERVATIONS OF CONVECTIVELY DRIVEN VORTEX FLOWS IN THE SUN 
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ABSTRACT 

We characterize the observational properties of the convectively driven vortex flows recently dis- 
covered on the quiet Sun, using magnetograms, Dopplergrams and images obtained with the 1-m 
balloon-borne SUNRISE telescope . By visual inspection of time series, we find some 3.1 x 1CV 3 vor- 
tices Mm min , which is a factor of ~1.7 larger than previous estimates. The mean duration of 
the individual events turns out to be 7.9 min, with a standard deviation of 3.2 min. In addition, we 
find several events appearing at the same locations along the duration of the time series (31.6 min). 
Such recurrent vortices show up in the proper motion flow field map averaged over the time series. 
The typical vertical vorticities are < 6x 10 -3 sec -1 , which corresponds to a period of rotation of some 
35 min. The vortices show a preferred counterclockwise sense of rotation, which we conjecture may 
have to do with the preferred vorticity impinged by the solar differential rotation. 
Subject headings: convection — Sun: granulation — Sun: photosphere — Sun: surface magnetism 



1. INTRODUCTION 

Solar surface convection is driven by localized down- 
drafts that collect the cold plasma returnin g to the solar 
interior after releasing internal energy (e . g ISpruit et al.l 
[19901: IStein fc NordlundlHoM Ir1^]ll998I) . Angular mo- 
mentum conservation forces the plasma to spin up as 
it approaches the sinkhole, and vortices are formed 
at the downdrafts. Such convectively driven vortices 
were theoret i cally p redicted and sought for long (e.g., 
ISpruit et all 119901: Ivan Ballegooijen et al.l Il998li. but 
their observational discovery is fairly recent dBonet et all 



20081 iWedemever-Bohm fc Rouppe van derVoortl 120091: 
Balmaceda et alj|2010t iGoode et al.ll2010l with the well- 
known earlier dete ction of a single large whirlpool by 
iBrandt et alJ lT988). In addition to supporting numeri- 
cal models of solar surface convection, the photospheric 
vortices may be of importance as heating sources for the 
outer solar atmosphere. The downdrafts not only advect 
vorticity but also magnetic fields, which arc intensified to 
kG field strengths in and around them. Buoyancy and 
the vertical geometry of the downdraft tend to align the 
magnetic field lines with the vertical, so that the spinning 
motions at photospheric levels can be propagated up- 
ward using the field lines as guides fe.g..lChoudhur ret al.l 
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119931: iZirkerl 119931 : Ivan Ballegooiien et alJll998j ). Waves 
thus excited transport photospheric energy that can be 
deposited in higher layers of the atmosphere. Moreover, 
downdrafts often trap structures of mixed polarity, so 
that the swirling motions wind up opposite polarity field 
lines, facilitating magnetic reconnection and the ensuing 
energy release. 

These convectively driven vortex flows are a recently 
discovered phenomenon poorly characterized from an ob- 
servational point of view. So far, we only know that the 
vortices are quite common, have no preferred sense of 
rotation at the sol ar equator, and last (at least) minutes 
(jBonet et al.ll2008h . They are also visible in the chromo- 
sphere, whe re they seem to be associated with significant 
blues hifts (|Wedemever-Bohm fc Rouppe van der Voortl 
1 2009ft. Most o f them are small-scale (< 0.5 Mm 
iBonet et al.ll2008f) , but so me have a much larger radius of 
influe nce (up to 20 Mm, lAttie et al.l 120091 : IBrandt et al.l 
119881 ). Lifetimes can be longer than 20 min, and several 
observables (such as circular polarization and G-band 
intensity) si multaneously indicate t he presence of vorti- 
cal motions ()Balmaceda et aLll2010[ ). In terms of global 
properties rather than individual eddies, the vertical vor- 
ticity inferred from proper motions seems to be higher 
in downflow re gions, suggesting excess vorticity in inter- 
granular lanes (Wane e t alJll995t iPotzi fc Braridll2q05| ) . 
iNisenson et al.l (|2003t ) searched for evidence of vorticity 
in the motions of isolated G-band bright points. 

Sunri se is a 1 -m balloon-borne solar tele- 
scope (|Barthol et al.l 120101 ) which, together with 
its Imaging Mag ne togra ph experiment (IMaX, 
iMartmez Pillet et al.l I2010T ) . provides time series 
with state-of-the-art high spatial resolution images and 
magnetograms. The dataset is ideal for a systematic 
characterization of the poorly known physical properties 
of the vortices. Thus this Letter presents a comprehen- 
sive observational characterization of the vortices in the 
photosphere. The actual data set and the procedure to 
detect vortices are described in § [5J The main results are 
summarized in § [3] Based on such results, we compare 
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the observed properties with the predictions of the 
numerical simulations of solar surface convection ( § [4j) . 

2. OBSERVATION AND ANALYSIS PROCEDURE 

The data were gathered with IMaX near the solar disk 
center on June 9, 2009 (UT 01:31-02:02; although the ex- 
act locati on is not known) dur i ng the first science fl ight of 
Sunrise (|Barthol et al.l[20Tol : ISofanki et al.ll20Toh . The 
IMaX magnetograph uses a LiNbOa etalon operating in 
double pass, liquid crystal variable retarders as the polar- 
ization modulator, and a beam splitter as the polariza- 
tion analyzer. We use here data recorded in the so - called 
V5-6 observing mode (see iMartmez Pillet et al.l 12010( 1 
where images of the four Stokes parameters were taken 
at five wavelengths along the profile of the magnetic- 
sensitive line Fe IA5250.2 A (±80, ±40 ml from line 
center, plus continuum at ±227 mA). After the science 
observing run a calibration set consisting of 30 in-focus 
and out-of-focus image-pairs was recorded for post-facto 
retrieval of the point spread function (PSF) us ing phase 
diversity (|Gonsalveslll982l : lPaxman et al.lll996f >. The sci- 
ence images were reconstructed by deconvolution using 
a modified Wiener filter and the calibrated PSF of the 
optical system. IMaX provides 85 mA spectral resolu- 
tion and between 0'.'15 and 0'.'18 angular resolution in 
the reconstructed images. Dopplergrams and magne- 
tograms are derived from th e Stokes parameters by usin g 
the approach described in ([Martinez Pillet et al.l[201Q] ). 
All in all, the reduction procedure renders time series 
of images, magnetograms, and Dopplergrams with a ca- 
dence of 33 sec, a spatial sampling of 0'.'055, and an ef- 
fective field-of-view (FOV) of 45" x 45". As inferred from 
the standard deviation of the polarization signals at the 
continuum wavelength, the circular polarization noise is 
5 x 10~ 4 in units of the continuum intensity. The ob- 
serving material analyzed here consists of a time series 
lasting 31.6 min. Movies were genera ted after rigid alig n- 
ment and jp-m o de su bsonic filtering ((Title et al.l ll989). 

iBonet et al.l (|2008| ) found the small-scale vortex flows 
by visual feature tracking of magnetic G-band bright 
points (BPs). The same downdrafts producing vortices 
also advect and concentrate magnetic flux (see § [1]), 
which often appears as BPs w hen the field strength 
is in the kG regime (see, e.g., iSanchez Almeida et al.l 
2004, and references therein). As Bonet and col- 
leagues acknowledge, the technique is rather limited since 
whirlpools without BPs are expected, and they escape 
from detection. Taking advantage of the combined high 
spatial resolution and high polarimetric sensitivity of 
IMaX/SuNRiSE data, we tried to detect and study vor- 
tices in longitudinal magnetograms, which are sensitive 
not only to kGs fields, but to plasmas with the full range 
of field strengths. In addition, we broaden the study 
using continuum intensity, line minimum intensity, line- 
of-sight (LOS) velocity, and line width (the last three 
parameters obtained from a Gaussian fit to the five sam- 
pled wavelengths). Individual vortices are identified and 
characterized through the following steps: 

(1) The detection is based on a visual inspection of 
longitudinal-magnetogram movies. Playing back and 
forth these movies, we identify those locations and time 
intervals where structures seem to rotate. 

(2) Once a vortex candidate is thus located, it is iso- 



lated in a 5"5 x 5"5 sub-field, where the corresponding 
sub-fields in the other four physical parameters are visu- 
ally inspected for swirling motions (see Fig. [IJ. 

(3) Horizontal velocity maps of the event in all the 
five parameters are created from proper motions. The 
horizontal motions are measured in these reduced sub- 
fields empl oying the local correl a tion t racking (LCT) al- 
gor ithm of iNovember fe Simonl (|1988l ). as implemented 
by iMolownv-Horas fe Yil (|1994D . and with a Gaussian 
tracking window of about 0"4 FWHM. In order to help 
the algorithm, the original images are interpolated in 
time and space so as to have a pixel of 0"028 and a ca- 
dence of Usee. The horizontal velocities obtained by 
comparing successive images are time averaged the du- 
ration of the event. Examples of such velocity maps are 
shown in Fig. [2j The size of the tracking window was 
chosen as a trade off to be large enough for the LCT al- 
gorithm to have structures to track, yet small enough to 
minimize the presence of several structures with differ- 
ent velocities. If the velocity maps do not show a regular 
closed shape in at least two physical parameters, then we 
discard the vortex candidate and start from step 1. 

(4) We compute the vertical vorticity, the divergence, 
and the curvature corresponding to the LCT horizontal 
velocities. Given the velocity U, the vertical vorticity 
(V x U) z can be interpreted in terms of the local angular 
velocity since a plasma in pure rotational motion has 
(V x U) 2 = 2ui, with w the angular velocity. Similarly, 
the curvature of such motion is k = 3 (V x U) Z |U| _1 , 
with At -1 the radius of curvature. Examples of vorticity 
and curvature maps are shown in Figs. [3j |4l 

(5) Using the LCT horizontal velocities, we track the 
evolution of passively a dvecte d tracers (corks) spread out 
all over the sub- field ( Yi 1992). As time goes by, the corks 
end up in the sinkhole, revealing the position of the sink 
according to different physical parameters (see the white 
points in Fig. [3]). If the sinkhole positions inferred from 
the different physical parameters do not agree within 0"4, 
then the vortex is discarded and we return to step 1. 

We find the curvature maps to be an efficient comple- 
mentary tool for vortex center detection. In most cases 
these maps show up the sinkholes as conspicuous point- 
like features (see Fig. [3J with positions that agree well 
with the centers determined as the final destination of 
the corks in the corresponding movie (c.f. Figs. OH]). 

The above list outlines the general procedure, but we 
do not disregard casual detections, e.g., when a second 
vortex was observed in any of the sub-fields corresponding 
to another vortex. Thus we find by chance some vortices 
which do not show up in the magnetogram signals. In 
addition to the LCT velocity maps of the individual vor- 
tices, we also computed the flow field for the full FOV 
during the full duration of the sequence. Isolated point- 
like features in the corresponding curvature map suggest 
the presence of vortices, and their existence motivates 
further inspection of magnetograms for swirling motions. 
As one can image from this cumbersome procedure, the 
FOV has been unequally searched. We focused on those 
regions were the magnetograph signals were largest, so 
that the effective FOV of our research is only 28"5 x 
28"5. This area is used in the estimates below unless 
otherwise stated. 

Two animations, one showing the event in Figs. [T0 
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and another with a different one, are given in the on-line 
material in the electronic edition of the Journal. 

3. RESULTS 

Following the procedure outlined above, we detected 42 
vortices with proper velocity maps. They imply a space- 
time-density of d ~ 3.1 x 10~ 3 vortices Mm -2 min -1 . 
In addition, 31 structures showing vortical motion were 
discarded because they did not fullfill our strict selection 
criteria. The selected events are used here to characterize 
the observational properties of vortices. 

We have assigned a duration to each vortex, computed 
as the time interval in which the vortex motions are clear- 
est. These durations span from 5min to 20min with a 
mean of r ~ 7.9 min and a standard deviation of 3.2 min. 
The interpretation of these intervals as life-times is not 
devoid of uncertainty, though. Often we shorten the in- 
terval to assure a most pure swirling motion. In addition, 
some vortices appear in the mean flow field correspond- 
ing to the full time series, indicating that they probably 
last longer than the time span of the series itself. These 
long lasting vortices often involve a complex behavior: 
several short-lived vortices appear and disappear in the 
same location, giving rise to recurrent vortices. The po- 
sition of their vortex centers may be static or drift with 
time. The recurrent vortices may or may not keep the 
same sense of rotation. In the latter case, however, we 
cannot consider the vortices to be strictly recurrent. We 
even find cases where the presence of a vortex is hinted 
as a clear pointlikc feature in the curvature image, but 
we failed to identify any vortex at that position during 
the sequence. 

The largest surprise of our analysis is the finding of 
a preferred sense of rotation: 27 counterclockwise vs. 
15 clockwise. Th is is a big difference with respect to 
iBonet et al.l ()2008l ) , where the two senses of rotation were 
observed equally. We analyze this issue in § |U 

Figure [3] shows histograms of vertical vorticities to 
characterize our measurements. The solid line corre- 
sponds to the vorticity in a region 2 Mm wide around 
a number of well defined vortices. The histogram con- 
sidering only 0.5 Mm is shown as the dotted-dashed line. 
These histograms reveal the signature of the vortices, 
which turn out to have vorticities up to 0.006 sec -1 , cor- 
responding to a period of rotation of some 35 min. These 
are vorticities inferred from the magnetograms, which 
are systematically smaller than those obtained from the 
other observables. The difference can be pinned down 
to the proper motion velocity field, which tends to zero 
outside the large magnetic concentrations, where the po- 
larization signals are low. The proper motions inferred 
from the other parameters extend throughout (Fig. [5]) , so 
that the histograms of vorticities have less contribution 
at low values and show extended tails (cf. the dashed 
line and the dotted line in Fig. [5J which represent the 
histograms of vorticities for the full FOV inferred from 
the magnetogram and the continuum intensity, respec- 
tively). 

If a vortex has its axis tilted with respect to the LOS, it 
should produce a characteristic Doppler signal similar to 
the rotation curve of a galaxy, with a close pair redshift- 
blueshift centered at the sinkhole. The expected signals 
are of the order of a few hundred ms -1 for moderate- 
high inclinations (30°), which are at the limit of our ob- 



servation. We unsuccessfully seek for such signals in the 
Doppler maps, meaning that the vortex motions are not 
highly tilted with respect to the horizontal plane. More- 
over, we note the discovery of hor izontal vortex flows n ear 
the edges of granules reported bv lSteiner et al.l (|2010D us- 
ing the same SUNRISE/IMaX data and that are likely to 
be of a different nature to those analyzed here. 

4. DISCUSSION 

Small-scale vortex flows in quiet Sun are detected using 
five different observational parameters: magnetogram, 
continuum intensity, line core intensity, LOS velocity, 
and line width. The fact that in most cases the detec- 
tion is consistent in three or more of these observables 
(showing different morphology) reinforces the reliability 
of the events found. 

The number den sity of vortices is ~1.7 times larger 
than that found by IBonet et al.l ()2008D . even though we 
have been far more strict here. The increase can be 
ascribed to the use of a larger variety of physical pa- 
rameters to detect the swirls. Most of the vortices are 
shortlived events observed during less than 10 min, but 
some of them last longer than the full time series, with 
recurrent vortices appearing in roughly the same place. 

Vortices have a typical vorticity smaller than 
6xl0 -3 sec -1 , which corresponds to a period of rotation 
of some 35 min. For ref erence, the large maelstrom found 
bv lBrandt et al.l (fl98 8) had a vorticity ten times smaller, 
with an associated period of some 6 hours. The mea- 
sured vorticities are generally much smaller than those 
predicte d by the numerical simu lations of magneto con- 
vection (jStein fc Nordlundlll998l Stein 2010 private com- 
munication). We think that the bulk of this difference 
can be attributed to the limited spatial-temporal resolu- 
tion of the observations. We are unable to identify vor- 
tices smaller than the tracking window, and/or lasting 
less than 8-10 frames, which sets an upper limit to the 
vorticity of some 0.04 sec -1 . Simulations indicate that 
the highest vorticities occur at the smallest resolvable 
scales and, thus, the predicted distribution critically de- 
pends on the numerical resolution of the simulation (see 
Fig. 31 in lStein fc NordlundllT998l ). 

The curvature maps (i.e., maps of the inverse radius 
of curvature) show the presence of vortices much better 
than the vorticity maps (see Figs. [3j H|). The vorticity 
is sensitive to the flow speeds, which are large outside 
vortices, creating spurious vorticity signals. The curva- 
ture, however, only enhances areas where the swirling 
motions occur at small scales, independently of the ve- 
locity. The sinkholes show up conspicuously as local ex- 
tremes in the curvature maps, and this new property 
should be employed when devising automatic algorithms 
to detect vortices. 

We find a preferred sense of rotation for the vortices 
(27 counterclockwise vs 15 clockwise). If the two senses 
of rotation were equally-probable then our observation 
would be highly unlikely (the probability is 4.4% assum- 
ing a binomial distribution). Nevertheless, the statistics 
is not large enough to provide a firm conclusion. The 
role of Coriolis forces on setting up this difference can 
also be discarded since the vortex motions involve time 
scales much too short to be affected by the solar rota- 
tion. The preferred sense of rotation may have to do 
with the solar differential rotation. The plasma pole- 
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Fig. 1. — A particular vortex as reflected in different observational parameters - from left to right: magnetogram, continuum intensity, 
line core intensity, LOS velocity, and line width. The panels show the average along the life-time of the event (~ 6.7 min), for every 
parameter. Horizontal scales are given in Mm. The 1 -Mm radius circles centered in the sinkhole are included for reference. 
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Fig. 2. — Horizontal velocity maps derived by LCT, from the proper motions of the parameters shown in Fig.l. The velocities are 
averaged over the life-time of the event (~ 6.7 min). Horizontal scales are given in Mm. The length of the black bar at coordinates (0,0) 
corresponds to 1.8 kms -1 . 

The 1 -Mm radius circles centered in the sinkhole are included for reference. 




Fig. 3. — Maps of vertical vorticity corresponding to the LCT horizontal velocities for the parameters shown in Fig.l. The white points 
indicate the final position of the corks in the cork movies (see text). Horizontal scales are given in Mm. The 1 -Mm radius circles centered 
in the sinkhole are included for reference. Vorticities are represented in the range ±6 X 10 — 3 s _1 , using a common grayscale in all the 
panels. 



ward from the sink tends to lag behind, whereas the 
plasma equatorward from the sink moves forward. Such 
a difference impinges a preferred counterclockwise sense 
of rotation in the northern hemisphere, and a clockwise 
sense in the southern hemisphere. Back-of-the-envelope 
estimates indicate that the effect produces the right or- 
der of magnitude vorticitvFl If this conjecture turns out 

10 Assume plasma separated by 40 Mm (i.e., the size of a su- 
pcrgranule) and converging to a sinkhole in the middle. Then the 
difference of differential rotations at latitude 30 deg produces a 
difference of velocities of some 35 msec -1 . If the circulation is 



to be corr ect, it natur a lly ex plains the difference with 
respect to iBonet et al.1 <|2008f >. whose observations cor- 
respond to the solar equator where there is no preferred 
sense of rotation. In our case the dominant counterclock- 
wise rotation is consistent with an observed FOV in the 
northern hemisphere. 

Thanks are due to R. Stein for discussions on the 

approximately conserved during the convergence, a vortex 0.5 Mm 
wide would have a vorticity of 0.01 sec -1 . 
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Fig. 4. — Curvature maps derived from the LCT horizontal velocities for the parameters shown in Fig.l. Horizontal scales are given in 
Mm. The 1 -Mm radius circles centered in the sinkhole are included for reference. Note how the curvatures clearly show the position of 
the sinkhole as an intense point-like feature. Curvatures are represented in the range ±5 X 10 -2 km^ 1 , using a common grayscale in all 
the panels. 
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Fig. 5. — Probability density function (PDF) of vertical vortic- 
ities obtained from LCT proper motions. The 5 min time aver- 
aged distributions for the full FOV are shown as a dotted line, 
and a dashed line, depending on whether they were derived from 
the continuum intensity or the magnetogram, respectively. The 
dotted-dashcd line corresponds to the vorticity in a region 0.5 Mm 
wide around a number of well defined vortices (with the vorticity 
signed so that all vortices have positive vorticity). The solid line 
also represents a local histogram considering a region of 2 Mm. The 
excess of vorticity at some 0.005 sec -1 is produced by the vortices. 
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